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APPLICATION  OF  THE  HERTZ  CONTACT  LAW  TO 
PROBLEMS  OF  IMPACT  IN  PLATES 


Prepared  by: 

Jackson  C.  S.  Yang  and  Do  Sup  Chun* 


ABSTRACT:  The  purpose  of  this  paper  is  to  determine:  (1)  the 
validity  of  t^e  static  treatment  of  pressure  on  thin  plates  for  the 
dynamic  cases  of  intermediate  impact  velocities;  (2)  the  use  of  the 
Hertz  law  for  the  contact  of  an  elastic  projectile  (steel)  on 
elastic  (steel),  viscoelastic  (Adiprene),  and  plastic  (Lexan)  bodies; 
(3)  the  use  of  viscoelastic  and  plastic  materials  as  shock  mitigators. 

Equations  for  stress  distribution  in  the  simply  supported  square 
plate  under  transverse  impact  have  been  developed  by  using  the  static 
theory  of  thin  plates  combined  with  the  Hertz  contact  law. 

Experimental  tests  were  performed  using  an  air  gun  to  impact  steel 
spherical  projectiles  on  steel  plates  with  and  without  Adiprene  and 
Lexan  coverings.  Two  spherical  projectiles  of  different  radii  were 
used.  The  radius  of  contact  between  the  projectile  and  the  covering 
was  measured  after  each  impact  test.  Strain  measurements  were  also 
taken  from  strain  gages  mounted  on  the  free  side  of  the  plate 
directly  under  the  point  of  impact.  Reasonable  agreement  between  the 
results  of  the  calculations  based  on  the  theoretical  analysis  and  the 
experimental  data  was  obtained. 
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two  bodies 
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INTRODUCTIOM 

During  the  past  decade,  increasing  attention  has  been  focused 
on  the  probless  attendant  to  the  collision  of  a  projectile  and  a 
target.  Thus,  the  Impact  problems,  as  well  as  the  fatigue  problems, 
have  been  investigated  and  developed  bv  many  authors  for  various 
conditions.  In  the  impact  studies,  the  Hertz  theory  has  been 
applied  even  though  some  questions  have  been  raised  concerning  the 
validity  of  the  theory  under  dynamic  conditions,  since  the  Hertz 
theory  is  based  on  the  quasi-static  contact  of  indenter  against  the 
infinitely  large  target.  Introducing  his  law,  Hertz  published  his 
original  analysis  of  elastic  contact  stress  under  static  loading  in 
1881  (Ref.  1),  and  his  analysis  revealed  that  the  Intensity  of 
pressure  q  over  the  surface  of  contact  between  two  general  curved 
surfaces  is  represented  by  the  ordinates  of  a  semi-ellipsoid  con¬ 
structed  on  the  surface  of  contact,  as  in  Figure  1;  thus 


in  which  q^  is  maximum  pressure  and  expressed  as 


q 


o 


3P 

Tv  a-*  b"" 


(2) 


where  P  is  total  applied  load,  and  his  force-indentation  law  is 


P 


■  n  a 


3/2 


(3) 


where  n  is  determined  with  the  properties  and  dimensions  of  two 
bodies,  and  a  is  the  relative  indentation. 

His  static  analysis  has  also  been  used  in  investigating  Impact 
of  elastic  bodies. 

As  an  example  of  the  impact  of  two  spheres  (Fig.  2),  if  Vj^ 
and  ^2  velocities  of  two  bodies  after  contact,  then 

4  -  Vj  +  Vg  (4) 
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With  the  above  relation,  using  equations  of  motion  of  two  bodies  and 
the  Hertz  force-indentation  law  (Eq.  (3)),  we  can  find  the  maximum 
compressive  force  P  acting  between  spheres  during  impact,  and  the 
corresponding  radius  of  surface  of  contact. 

Reference  S  indicates  that  for  the  case  of  impact  of  elastic 
bodies  at  model t.te  velocities,  the  problems  of  elastic  contact  and 
elastic  impact  are  in  essence  identical  (Love  (1934)),  and 
reference  3  indicates  that  the  quasl-etatic  treatment  is  found  to 
Approximate  very  closely  the  values  of  stress  found  in  dynamic 
impacts . 

C.  V.  Raman  (Ref.  4)  applied  the  Hertz  theory  (Eq.  (3))  to  the 
investigation  of  the  coefficient  of  restitution,  and  foimd  that  for 
moderate  thicknesses  of  plate  (0.138”  1”),  tb9  theoretically  cal¬ 

culated  and  experimentally  observed  values  agreed  well. 

In  further  investigations,  Werner  Goldsmith  (Ref.  5)  extended 
his  experimental  observations  to  intermediate  velocities  (30-300  fps) 
of  impact,  using  the  Hopklnson  bar  technique,  and  vindicated  the  use 
of  the  Hertz  law  for  hard  tool  steel  even  when  some  permanent  set 
was  produced;  however,  he  denied  its  application  to  aluminum. 

For  the  purpose  of  determining  the  validity  of  applying  the 
static  treatment  of  plates  and  Hertz's  law  to  thin  plates  under 
intermediate  velocities  of  impact,  the  present  report  gives  an 
analytical  solution  for  the  stress  distribution  of  the  plate  under 
transverse  impact.  This  is  done  by  the  combined  application  of  the 
static  theory  of  plates,  together  with  Hertz's  law. 

A  comparison  is  then  made  with  the  experimental  data  obtained 
using  an  air  gun  for  intermediate  velocities  of  impact.  Two  pro¬ 
jectiles  with  different  radii  of  curvature  were  used  to  impact 
against  steel  plates  of  various  thicknesses  with  and  without  poly- 
urethene  (Adlprene)  and  polycarbonate  (Lexan)  coverings.  The  shock 
mitigation  properties  of  the  two  claddings,  Adlprene  and  Lexan,  were 
also  obtained. 

THEORETICAL  ANALYSIS 


For  the  case  of  impact  between  two  spheres  of  mass  m.  and  m^, 
investigations  show  that  the  duration  of  impact,  l.e.,  the  time 
during  which  the  spheres  remain  in  contact,  is  very  long  in 
comparison  with  the  period  of  lowest  mode  of  vibration  of  the  spheres 
(Ref.  6).  Vibrations  can  therefore  be  neglected,  and  it  can  be 
assumed  that  the  force-displacement  relation  established  for  static 
conditions  holds  during  impact.  The  force  is  basically  a  power  law 
given  by 


P 


n  a 


3/2 


(5) 
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where  P  is  the  compressive  force,  and  the  constant  n  Is  given  by 


n 


4 _ 

3»  (K^  +  Kg) 


(6) 


Rg  -  Radii  of  the  spherical  surfaces  of  the  two  bodies  at 
the  point  of  contact 


K 


1 


1  - 


2 


w  E 


2 


in  which  E  and  v  are  Young's  modulus  and  Poisson's  ratio, 
respectively. 


This  is  similar  to  examining  the  material  as  a  nonlinear 
spring.  In  our  problem,  a  spherical-nosed  projectile  is  Impacted 
against  a  rigid  flat  surface.  Then  equation  (6)  may  be  simplified 
to 


n 


4  R, 


1/2 


3w(Kj  +  Kg) 


(7) 


in  which  subscripts  1  and  2  represent  plate  and  projectile, 
respectively. 

During  the  impact  of  such  a  system,  the  following  differential 
equation  may  be  written: 


•*  3/2 

fflg  a  +  n  a  -  0 


(8) 


The  solution  to  the  system  differential  equation  gives  the 
displacement,  velocity,  and  acceleration  of  the  spherical-nosed 
projectile  with  respect  to  time.  With  the  displacement  of  the 
spherical  projectile  known,  the  force  acting  on  the  plate  at  any 
time  can  be  obtained  from  equation  (5). 
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Plate  of  Infinite  Thickness 


The  maximujD  compression  of  the  cladding  and  the  duration  of 
impact  can  be  determined  directly  by  requiring  the  work  done  by 
the  cladding  to  equal  the  kinetic  energy  of  the  system.  The 
deflection  of  the  plate  is  neglected,  since  it  is  very  small  whan 
compared  to  the  deflection  of  the  cladding.  In  this  case  (Ref.  6) 


Hertz's  contact  law  is  also  used  to  estimate  the  maximum 
bearing  stress  which  the  system  must  sustain.  The  stress  distri¬ 
bution  is  hemispherical  in  shape  over  a  circle  of  contact  of 
radius  c  (Ref.  6)  given  by 


r  3  w  P  (Kj  +  Kg)  RjRg  1 
L  4  (Rj  +  Rg)  i 


1/3 


For  Rj^  •  • 


-|’3/4  w  P  (Kj  +  Kg)  Rgj 


1/3 


We  have  the  maximum  compressive  force  as 


p  .  „  c.  3/2 . 

max  max 


w  mg 


3/5 


3  w  (Kj  +  Kg) 


2/5 


(11) 


(12) 


(13) 


and  the  corresponding  radius  of  the  surface  of  contact  can  be 
obtained  by  substituting  equation  (13)  into  equation  (11) , 


c 


max 


)2/5  ^  15t  (K^  Kg)  mg  v  ^ 


(14) 
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If  the  two  bodies  are  of  the  same  material, 

3/5 


4(R2) 


1/5 


max 


(tI  *  *2''^) 

n  '2/5' 


3v 


'2(1  -  v^)  ' 
.  w  E 


(15) 


S»ax  “  ^**2' 


2/5 


[ 


15  (1  -  V 


8  E 


)  °>2  ^ 


2  nl/5 


(16) 


The  peak  bearing  stress  is  obtained  by  equating  the  sum  of 
the  pressures  over  the  contact  area  to  the  compressive  force  P. 
For  the  hemispherical  pressure  distribution  this  gives 


3  P 


max 


2  r  c^ 


max 


(17) 


This  stress  can  be  written  as 


a  function  of  displacement 


^o 


E 

t(l  -  v^) 


(18) 


With  the  maximum  compression  force,  And  the  maximum  radius 

of  the  surface  contact,  c  ,  as  equations  (13)  and  (14)  for 
different  materials,  and  equations  (15)  and  (16)  for  the  same 
materials  known,  the  stress  can  be  calculated  for  the  outer  sur¬ 
face  of  the  plate  at  the  pole  by  the  thin  plate  theory  (Ref.  7). 

Plate  of  Finite  Thickness 


In  the  previous  analysis,  due  to  the  cladding  the  deflection 
of  the  plate  is  neglected  since  it  is  very  small  when  compared  to 
the  deflection  of  the  cladding.  However,  for  impact  tests  without 
cladding,  the  bending  of  the  plate  should  be  taken  into  considera¬ 
tion. 


From  reference  7  by  Goldsmith,  a  differential  equation  of 
motion  was  written  for  an  approximate  solution  to  the  case  of  the 
central  transverse  impact  of  a  sphere  on  a  rectangular  plate  simply 
supported  along  all  edges. 
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.2  , 

^  i  F(a)  + 


L6  L 


3(1  -  v^) 
P  K 


F  (a)  -  0 


The  elastic  force-indentation  lav,  equation  (5),  is  now 
enployed  to  eliminate  the  term  F(a)  from  the  equation.  Intro¬ 
ducing  dimensionless  variable  defined  by 

5  .  -2_  .  I  .i 


equation  (9)  may  be  written  as 


vitb 


a  -  0  , 


dt  J 


3/2  ^ 


1  for  t  -  0 


where  f,  a  constant  with  dimensions  of  time,  is  so  chosen  that  the 
first  two  coefficients  in  the  transformed  equations  are  unity.  For 
a  striking  sphere  of  radius  R,  the  parameter  \  is  given  by 


pjd  - 


1  -  V, 


®2 

"TT  +  . . 


1  .  Vi  I  -  Vj 


where  subscripts  1  and  2  denote  the  plate  and  sphere,  respectively. 

The  results  of  numerical  integrations  of  equation  (21)  for 
several  values  of  1  are  shown  in  Figure  3  for  dimensionless  dis¬ 
placement  histories  as  a  function  of  dimensionless  tlme/-u 


From  Hertz's  contact  force  law 


'(»vj 


this  gives  us  the 
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force  F  acting  on  the  plate  at  anytime »  having  taken  into  considera¬ 
tion  the  bending  of  the  plate. 

Derivation  of  Stress 


1.  Case  of  hemispherically  distributed  contact  pressure. 

According  to  Hertz,  the  intensity  of  pressure  q  over  the 
surface  of  contact  is  represented  by  the  ordinates  of  a  hemisphere 
of  radius  c  constructed  on  the  surface  of  contact  as  shown  in 
Figure  4  (special  case  of  Fig.  l)(Refs.  6,  8). 

Taking  a  coordinate  system  as  shown  in  Figure  4,  the 
deflection  of  the  plate  must  satisfy  the  differential  equation, 


d^w  o  3^w 

— I  ^  — 9 — 2  — I 


where  D  - 


Eh* 


12  (1  -  v^) 


q(x,y) 

D  ■ 


(23) 


(flexural  rigidity  of  the  plate) 


q  (x,y) 


(y  -  n)^ 

- 2 - 


Solving  the  differential  equation  with  the  appropriate 
boundary  conditions,  we  have  the  general  form  of  the  deflection  of 
a  simply  supported  thin  plate  subjected  to  the  transverse  impact 
load  as  follows  (Hef.  7): 


W  - 


Amn 


^ 


sin 


mvx 

a 


sin 


nvy 

IT 


(24) 


where  A 


mn 


GF 


abvc‘ 


rn 


-  r 


sin 


rav(5+ycose) 


sin  nT(n-fysln6) 
- B - 


rdrd6  (25) 


By  the  assumption  of  a  thtr.  plate,  o  «  0,  we  obtained  the  general 
forms  of  stresses  as  follows: 
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% 

% 


12z  p 
7T;Z  mSl 

tl  T 


-  12z(l  -  v)  p  P  Amn  m»x  cos  nzy 

V 


For  a  square  plate  (a  «  b  and§-q-  a/2), 


2  <D  CD  2  2 

12a  z  p  p  Aina(m  +vn  )  ra»x  nxy 

"*  ■  73-2-  m5l  n£l  TTT - STT  “a^ 

hr  (m  +  n  ) 


,  .  S,  L,  sin  !!2i  sin  HU 

y  .3  2  m*l  n“l  ,  3  2.2  a  a 

■'hr  (m  +  n  ) 


12a"z(l  -  v)  p 

“■3“^ -  m-1 

h  X 


Amn  (■)  (p) 


(m 


27? 

n  ) 


cos 


mxx 

a 


cos 


nxy 

a 


Since  maximum  bending  moments  occur  at  the  center  of  the 
(x  <-  y  and  the  maximum  bending  stresses  at  Z 


(26) 


(27) 


(28) 


(29) 


(30) 


(31) 

plate 
we  have 
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a  )  m  „ 

x'ms  .22 
n  T 


2  - 


2  2 

y;  ^  Amnlm  ♦  to  )  ^ 

*0,5  0-173,5  (.2'  "  ^7j2 


(-=4-2-  -  1) 


(32) 


®y^»«x 


6o*  "  "  ABn(vB  +  0  )  y  ,v 

.-0,5  0-0,5  “,2";  ;2)2 


(^4-^  -  1) 


(33) 


T  -  0 

xy 


(34) 


Notice  that  o  )  -  o  ) 

z'aaz  y'maz 


2.  Caee  of  uolforaly  distributed  coot act  pressure « 

Y.  M.  Tsai  aod  R.  Kolsky  (Ref.  3)  fouod  io  the  study  of 
wave  propagatloo  that  the  theoretical  aoalysis  gave  close  agreeMOt 
with  the  ezperlMotal  pulse  shapes,  by  assuaiog  that  the  area  of 
cootact  has  a  flolte  radius  aod  stress  is  distributed  uoiforaly 
over  this  circular  area.  Therefore,  for  the  siaplificatioo  of  the 
problea,  «e  cao  assuae  that  the  coapressive  force  is  uoiforaly 
distributed  over  the  cootact  area  as  shovo  io  Figure  5.  Sioce  the 
coapressive  force  P  is  uoiforaly  distributed  ovt^r  the  cootact  area 


Aao 


"jy  r.»tg9)  drde 


(35) 


Provided  that  the  circle  r  -  c  reaaios  eotirely  ioside  the  bouodary 
of  the  plate,  the  evaluatioo  of  the  iotegral  equatioo  (31)  gives 
the  express loo 


Aao  -  - J,  (Y  C)  sio  sio 

abc  Y.. 

mu 


(36) 


io  vhlch  Y, 
tic  of  order  ooe. 


aod  J. 


with  the  arguaeot  y 


(Y _ C)  is  the  Bessel  fuoc- 

ao 

(Ref.  11).  Therefore, 

■Q 
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the  expression  for  the  deflection  of  the  plate  becomes,  by  sub¬ 
stituting  equation  (36)  into  equation  (24) 


W  - 


8F 

“5 - 

T  Dabc 


m-J , 3 , 


5  n*lj^,5  r~7T?\~  /  ..^\  t  5/2 


sin  ^  sin  BIZ 
a 


(37) 


For  a  square  plate  a  -  b  and  ti  -  5  “  a/2 

3  OD  Ok 


W  - 


8Fa* 

»^Dc 


5  n-1.3.5  ”n^)^77- 


(=4-^  -1) 


sin  ^  sin 


(38) 


As  in  the  previous  section,  the  stress  at  the  center  of  the 
surface  of  the  plates 


a  )  “  o  ) 

x'max  y'^max 


48Fa 


ro*j, 3,5  n“*l^^,5 


r  /  «\  /  1  vra-fn  .  2  2. 

Jl(Y^C)(-l)  vni  -t^vn  ) 

- - T - 27572 - 

(m  -I-  n  ) 


(39) 


T  -  0 
xy 


EXPERIMEIITAL  TESTS 


Test  Apparatus 


An  air  gun  and  two  steel  projectiles  having  0.5-inch  and 
1.0-inch  radii  spherical  heads  were  used.  Air  pressures  up  to 
60  psi,  velocities  up  to  85  fps  were  obtained. 

Test  Specimens  and  Coverings 


Steel  square  plates  (6  in.  x  6  in.)  in  various  thicknesses, 
and  Adiprene  (h  -  0.29  in.)  and  Lexan  (h  -  0,19  in.)  coverings 
were  used.  Rockwell  hardness  tests  were  carried  out  to  obtain 
the  mechanical  properties  of  steel,  and  the  results  were  interpreted 
using  the  ASM  Handbook  (Ref.  10). 
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The  mechanical  properties  of  the  specimens  and  coverings  are 
as  in  Table  1. 


Comparison  of  the  Test  Results  and  the  Computed  Values 


The  comparison  of  experimental  results  with  theoretical  values 
for  velocity  of  Impact  versus  radius  of  contact  for  Lexan  and 
Adlprene  coverings  is  shown  In  Figures  8  and  9,  respectively. 
Experimental  and  theoretical  values  for  strain  versus  velocity  of 
Impact  on  the  steel  plate  with  and  without  coverings  are  shown  In 
Figure  10  for  1-lnch  radius  projectiles,  and  Figure  11  for  1/2-lnch 
radius  projectiles. 

DISCUSSION  AMD  CONCLUSION 


The  radius  of  contact  between  the  projectile  and  the  covering 
wao  measured  after  each  Impact  test.  These  values  were  compared 
with  the  theoretically  calculated  values  using  the  Hertz  law  for 
the  contact  of  elastic  (projectile)  and  viscoelastic  (Adlprene) 
bodies,  and  elastic  (projectile)  and  plastic  (Lexan)  bodies.  It 
was  found  that  there  was  a  reasonable  agreement  between  experimental 
and  theoretical  values.  This  confirms  the  validity  of  using  the 
Hertz  contact  law  In  problems  of  impact,  and  It  also  supports  the 
conclusion  that  the  application  of  the  Hertz  law  can  be  extended  to 
the  contact  of  viscoelastic  bodies  (Refs.  5  and  9). 

In  order  to  compare  the  experimental  results  with  the  theoretical 
values  of  the  strain  for  the  plate  with  coverings,  compressive  forces 
and  radii  of  contact  between  projectiles  and  coverings  were  first 
calculated.  Then  the  plate  was  considered  to  be  subjected  to  these 
forces  without  coverings.  The  values  of  stress  and  rtrain,  using 
the  biaxial  stress-strain  relations,  were  calculated  by  the  computer. 

The  experimental  data  were  obtained  from  the  strain  gages, 
mounted  on  the  free  side  of  the  plate  directly  under  the  point  of 
impact  and  connected  to  an  oscilloscope. 

From  the  comparison,  it  was  found  that  the  experimental  results 
were  In  good  agreement  with  the  theoretical  values  for  the  steel 
plate,  considering  the  fact  that  there  existed  a  certain  depth  of 
layer  of  covering  between  projectile  and  plate,  and  a  possible 
friction  loss  between  the  covering  and  the  plate. 

For  the  different  radii  of  projectiles,  it  was  found  that,  even 
with  the  same  mass  and  at  the  same  velocity  of  impact,  the  projectile 
of  larger  radius  produced  a  greater  compressive  force,  as  well  as  a 
larger  radius  of  contact  between  bodies,  thus  producing  a  larger 
strain,  both  theoretically  and  experimentally. 

From  an  evaluation  of  the  theoretical  results  for  both  the 
hemlspherlcally  and  uniformly  distributed  contact  pressure,  the 
latter  was  found  to  be  closer  to  the  experimental  results.  However, 
since  all  of  the  results  are  reasonably  close,  both  cases  can  be 
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used  to  predict  the  approximate  Impact  characteristics  of  thin 
plates. 

It  was  noticed  that  under  the  above  conditions  all  plates 
tested  remained  with  the  elastic  range. 

From  the  comparison  of  steel  plate  with  and  without coverings, 

It  was  found  that  Adlprene  served  as  the  best  mltlgator  of  shock. 

It  reduced  the  magnitude  by  a  factor  of  10  on  the  average.  Lexr.< 
also  served  as  a  good  mltlgator,  where  a  factor  of  5  was  obtalued. 

In  addition,  experimental  tests  were  also  carried  out  on  the 
plates  without  coverings  where  plastic  yleldlngs  (a  slight  permanent 
set)  occurred  at  the  point  of  contact.  Observing  the  results,  the 
same  conclusion  can  be  obtained  that  the  Hertz  law  Is  applicable  for 
steel  plates.  This  also  verifies  the  conclusion  reached  by 
Werner  Goldsmith  (Ref.  5)  that  the  Hertz  law  Is  valid  under  proper 
circumstances  even  when  some  permanent  set  Is  produced. 

It  Is  further  suggested  that  the  problem  of  large  plastic 
deformation  be  studied,  combining  the  theory  of  plasticity  with 
some  modification  of  the  Hertz  law. 
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FIG.  1  CONTACT  PRESSURE  DISTRIBUTION, 
TWO  GENERAL  CURVED  SURFACES 


FIG.  2  IMPACT  OF  TWO  SPHERES 
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FIG.  8  VELOCITY  OF  IMPACT  VS  RADIUS  OF  CONTACT 
(LEXAN  COVERING) 
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FIG.  10  STRAIN  VS  VELOCITY  OF  IMPACT  (  V'-R  PROJECT) 
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